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ABSTRACT 

We present optical light curves of twenty-nine novae in M31 during the 2009 and 2010 observing 
seasons of the Palomar Transient Factory (PTF). The dynamic and rapid cadences in PTF monitoring 
of M31, from one day to even ten minutes, provide excellent temporal coverage of nova light curves, 
enabling us to record the photometric evolution of M31 novae in unprecedented detail. We also detect 
eight of these novae in the near ultraviolet (UV) band with the Galaxy Evolution Explorer (GALEX) 
satellite. Novae M31N 2009-10b and 2010-lla show prominent UV emission peaking a few days prior to 
their optical maxima, possibly implying aspherical outbursts. Additionally, our blue-shifted spectrum 
of the recent outburst of PT And (M31N 2010-12a) indicates that it is a recurrent nova in M31 and 
not a dwarf nova in the Milky Way as was previously assumed. Finally, we systematically searched 
for novae in all confirmed globular clusters of M31 and found only M31N 2010-10f associated with 
Bol 126. The specific nova rate in the M31 globular cluster system is thus about one per year which 
is not enhanced relative to the rate outside the globular cluster system. 
Subject headings: surveys - nova, cataclysmic variables 



1. INTRODUCTION 

Classical novae (CNe) are interesting to study for at 
least two reasons: first, they provide a laboratory to 
investigate ther monuclear runaways in semidegenerate 
conditions (e.g., iBode fc 'Ev ans 2008); second, they have 
been proposed as prog enitors of type la supernovae (e.g., 
iStarrfield et al.1 fl988l ). Observations of Galactic novae 
are limited due to extinction and distance determina- 
tion. CNe are frequent but faint relative to supernovae. 
Hence, people pursue novae in the nearest galaxies. 

CNe in M31 (The Andromeda Galaxy) have been 
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(jDarnlev et al.|[200l . 

The Palomar Tran sient Factory (PTF) (|Law et al.l 
2009; Rau et al.ll2009f ) acquired nearly daily monitoring 
of M31 in the R band during the 2009 and 2010 observ- 
ing seasons. It provides three advantages for the study 
of novae in M31: first, the large field of view (FoV) of 
7.2 square degrees covers the main region of M31 in one 
snapshot and with a spatial resolution of l'.'O. Second, 
PTF has fast and dynamic cadences from one day to 
as short as ten minutes. This temporal coverage gives 
unprecedentedly detailed light curves. Third, the limit- 
ing magnitude is deeper than 20 mag (a little worse near 
the center of M31 where the background is strong and 
varies quickly). Henc e, given the mean di stance modu- 
lus of 24.36 for M31 (|Vilardell et al.|[2010h . unless some 
very faint novae are extincted heavily by local dust, PTF 
should see nearly all of the novae (weather-permitting) 
in M31 and follow them up down to a magnitude of 
M R 4.4. 

These data also afford the opportunity to explore 
known recurrent novae (RNe) in M31 and perhaps dis- 
cover new one s. According to the catalog maintained by 
IPietschl (|2010( ). about thirteen recurrent novae have been 
discovered in M31. RNe are possible progenitors of Type 
la supernovae, since their white dwarf mass is thought to 
be cl oser to the Chandrasek har mass limit than ordinary 
CNe (|Starrfield et al.|[l98l 

Another issue related to novae that can be assessed 
with PTF archives is whether the nova rate in the globu- 
lar cluster (GC) system is enhan ced relative to that out - 
side the globular cluster system. ICiardullo et a l. (1990b) 
did a search f or novae in fifty- f our M 31 G Cs based on the 
Ha sur vey of ICiardullo et al.l (|1987t ) and ICiardullo et aTl 
Ijl990al) . In a mean effective survey time of approximate 
two years, they did not find any nova outbursts in any of 
the M31 GCs. ITomanev et al.1 (|1992f) reported another 
one-year search of over two hundred M31 GCs with Ha 
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TABLE 1 

Observation summary 



PTF 



GALEX 



Filter 

^effect (A) 
AA (A) 

^expose (s) 
FoV (deg 2 ) 
Cadence 
Epoch 



Mould R 
6581 
1251 
60 
7.2 

lOmin to days 
726 



NUV 
2316 
1060 
740-1703 

1.2 
1-2 days 
51 (the center field) 
22-27 (outskirt fields) 



emission but did not find any nova outbursts, providing 
an upper limit on the rate of 0.005 nova yr _1 GC -1 . 
The only GC nova fou nd in M31 is M31N 2007-06b 
(jShafter fc Quimbvl[200l . iHenze et all (|2009eD then de- 
rived a rate of 1.1 nova yr _1 for the M31 GC system, 
which is comparable to the rate (per unit mass) outside 
the globulars. A systematic search for novae in M31 GCs 
using the PTF archives is presented here and also sup- 
ported this conclusion. 

In a nova outburst, except for a short "fireball" phase 
where radiation is given out by the thermal emission 
of the photosphere, UV photons are thought to orig- 
inate at the surface of the white dwarf (WD), deep 
within the nova envelope, while optical photons are 
released near the surface of the envelope and result 
from absorption and re-emission of relatively low ion- 
ization species such as Cr II, Fe II and other heavy met- 
als generated during the outburst. Nova shells there- 
fore serve as passive photon converters whose spec- 
tra reflect the reprocessing of incident UV light to 
longer wavebands. Eventually, the expansion decreases 
the temperature and dilutes a bsorbing species to fi- 
nally reveal the UV photons (iKato fe Hachisul 119941 : 
iGonzalez-Riestra k, Krautterl I1998D . Therefore, to ex- 
plore the UV light powering the nova emission, we 
searched in parallel for any contemporaneous ultra-violet 
observations of CNe in M31 with the GALEX satellite 
(jMartin et al.1 120051 ). Finally, we obtained near UV- 
optical light curves for eight novae in M31. 

This paper is arranged as follows. Section [5] describes 
PTF and GALEX observation and data reduction as 
well as data from other sources. Section [3] presents light 
curves of classical novae and describes their morphology. 
Section 0] discusses novae in globular clusters. Section [5] 
focuses on recurrent novae. Section [6] discusses ordinary 
and bizarre novae observed in UV band. Section [7] revis- 
its the MMRD relation with our novae sample. Finally 
we conclude this paper in Section [H 

2. OBSERVATION 

2.1. Nova Sample 

We examined all novae as reported bv lPietschl (pOl O'FI 
in the PTF image archives and obtained the light curves 
of twenty-nine novae. Among them, eight novae also 
have UV light curves from GALEX. Table Q] gives main 
parameters of PTF and GALEX observations. Table [2] 
describes the nova sample and Figure Q] illustrates their 
spatial distribution. 

11 http: / /www. mpe.mpg.de/~m31novae/opt/m31/index.php. 
This site maintains a comprehensive catalog of novae in M31. 



2.2. PTF 

PTF (jLaw et alJl2009h employs the 1.2-m Oschin Tele- 
scope at Palomar Observatory with an array of twelve 
CCD chips, each having 4Kx2K pixels. One of the chips 
is dead for unknown reasons. The pixel size is l'.'O/pixel 
and thus the total FoV is 7.2 square degree. 

2.2.1. Survey Observation 

The time coverage of PTF observations of M31 during 
the 2009 and 2010 seasons is shown in Figure [5] The 
first observing season is from 2009 August 18 to 2010 
January 17 and the second is from 2010 June 30 to 2010 
December 15. From 2010 August 28 to 2010 September 
7, PTF observed M31 with a fast cadence of about ten 
minutes. All PTF images were taken with a R-band fil- 
ter. The average seeing during the two seasons is 2'.'5 
with a median value of 2'.'4. We removed less than 1% of 
the images with seeing worse than 3'.'5, because the high 
stellar density in M31 leads to strong overlaps among 
adjacent objects under bad seeing conditions. Our final 
image database samples over seven hundred epochs on 
7987 individual CCD images. 

2.2.2. Data Reduction and Photometry 

All images were processed through the PTF LBL 
pipeline (see lLaw et aLl [2009), obtaining astrometric 
solutions of l'.'O accuracy. The calibration of relative 
photometry proceeded as follows. We first use SEx- 
tractor to make a source catalog for each image. We 
then choose the catalog with the largest number of 
14 mag < R < 17 mag sources as reference. Then, other 
non-reference catalogs are spatially matched to the ref- 
erence using source positions and forced to agree with 
the reference photometry after non-constant sources are 
removed with a standard outlier procedur43- The result- 
ing photometric uncertainty is usually around 0.05 mag. 
The absolute calibration is achieved by matching the 
reference sources to sources in the catalog presented in 
iMassev et all pOOl ). 

Photometry is also done with SExtractor. Background 
is estimated in a region of 20'.' x 20 7 .' centered at each tar- 
get. Fluxes are measured within a Kron-like elliptical 
aperture. The first order moment of e ach object deter- 
mines the Kron radius r (jKronl 119801 ) and the second 
order moment determines ellipticity and position angle. 
The aperture has an area of 6.257rr 2 . For a null detec- 
tion, a 3-<7 limiting magnitude is estimated. 

In order to assess the quality of photometric results 
from SExtractor, we constructed light curves of three 
novae (M31N 2009-09a, M31N 2009-10a and M31N2009- 
11a) with three methods: 1) the one we described 
above; 2) aperture photometry with an aperture radius 
of seeing and with constant local background; 3) aper- 
ture photometry with an aperture radius of seeing and 
with local background estimation by linear regression. 

12 NIST/ SMATCH e-Handbook of Statis tical Methods at 
http : / / www.itl . nist . gov / div898/handbook / 
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Fig , 1. — The spatial distribution of the nova sample in M31. The upper panel is an image of M31 taken by GALEX (Gil dc Paz et al. 
120071) . The bottom panel is a mosaic image of M31 from PTF (Credit: P. Nugent). The dashed box in the upper panel shows the physical 
size of the bottom one. The actual FoV is larger than the size of this image. The red circles denote positions of novae only observed in 
PTF. RNe are highlighted by the yellow circles. The magenta circles denote novae observed in both PTF and GALEX. 
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TABLE 2 

Sample of Novae 





Oi 







Discovery Date 


Alii^iiitLidc 




Spectral 




IV ova, Neiiiic 


(J2000) 




(J2000) 


by Others 


at Discovery 


Date 


TvDe 


Reference 


M31N 2009-08b 


00 ?l 44 m 09 1 ' 


87 


+41°48'50('8 


2009/08/09.780 


17.2(H) 1 


2009/08/18.423 




i 


M31N 2009-08e 


00 h 42 m 36 s 


23 


+41°18'01?6 


2009/08/25.896 


17.9(R) 


2009/08/26.400 


Fell 


2 


JV1311N 2009-09a 


( i ( i / ( id/fi lies 

00 42 2o 


08 


+41 04 01. 


2009/09/02.078 


17.1(U) 


or\nn /no /1 o a c a 

2009/08/18.454 


Fell 


3 


M31N 2009- 10a 


00 h 45 m 14 s 


01 


+42°04'38?8 


2009/10/03.619 


17.1(U) 


2009/09/30.494 


Fell 


1 


M31N 2009-10b 


00 h 42 m 20 s 


77 


+41°16'44'/5 


2009/10/11.414 


18.8(R) 


2009/10/10.376 


Fell 


5 


M31N 2009- 11a 


00 h 43 m 04 s 


76 


+41°4lW2 


2009/11/03.559 


17.8(U) 


2009/10/31.321 


Fell 


6 


M31N2009-llb 2 


00' l 42 m 39 3 


61 


+41°09'03('2 


2009/11/06.523 


18.6(R) 


2009/11/07.199 




7 


M311N 2009-lla 


c\C\h A A vn ics 
00 44 lb 


85 


+41 lo Oo. 


nnnn / i i /in lO/l 

2009/11/ 19.194 


in 1 /OA 

18.1(H) 


2009/ 12/01.201 


rcll 


8 


M31N 2009-lle 


00 ,l 42 m 35 s 


31 


+41°12'59'/1 


2009/11/21.138 


18.8(Swift uvw2) 


2009/12/01.201 


Fell 


9 


M31N 2009-12a 


00' l 40 m 19 s 


40 


+41°15'47?6 


2009/12/22.488 


16.6(U) 


2009/12/10.109 




10 


M31N2010-01a 3 


00' l 42 m 56 s 


74 


+41°17'21'/0 


2010/01/11.13 


17.6(R) 


2010/01/11.130 


Fell 


11 


M31N 2010-06a 


00' l 43 m 07 s 


56 


+41°19'49('0 


2010/06/28.014 


18.1(R) 


2010/06/30.397 


Fell 


12 


M31N 2010-Oob 


00 44 22 


46 


i A 1 O o a f 1 At t r 

+41 28 14. 5 


2010/06/28.014 


19.1(R) 


2010/06/30.353 


Fell 


13 


M31N 2010-06c 


00 h 44 m 04 s 


48 


+41°28'34?2 


2010/06/26.084 


17.8(R) 


2010/06/30.353 




14 


M31N 2010-06d 


00' , 42 m 55 s 


61 


+41°19'26'/0 


2010/06/24.02 


19.5(Swift uvwl) 


2010/07/04.397 


Fell 


15 


M31N2010-07a 


00 h 43 m 20 s 


11 


+41°21'23'.'7 


2010/07/07.036 


20.6(Swift uvwl) 


2010/07/07.328 


Fell 


16 


M31N2010-09a 


00 ,, 42 m 23 s 


32 


+42°17'08'.'6 


2010/09/04.552 


17.2(U) 


2010/09/01.372 


Fell 


17 


M31N 2010-09b 


00' , 43 m 45 s 


53 


+41°07'54'.'7 


2010/09/30.412 


17.7(R) 


2010/09/30.230 


Fell 


18 


M31N 2010-09c 


00' 1 38 m 09 s 


06 


+40° 37' 25'.' 9 


2010/09/03.30 


18.7(U) 


2010/08/13.235 


Fell 


19 


M31N2010-10a 


00 ,l 42 m 45 s 


84 


+41°24'22'.'2 


2010/10/05.551 


17.6(R) 


2010/10/08.163 


Fell 


20 


M31N 2010-10b 


00 ,t 42 m 41 s 


51 


+41°03'27'.'3 


2010/08/19.055 


18.9(R) 


2010/10/08.163 


Fell 


21 


M31N2010-10c 


00 ll 44 ra 26 s 


56 


+41°31'13'.'8 


2010/10/13.557 


17.8(U) 


2010/10/13.113 


Fell 


22 


M31N2010-10d 


00' , 42 m 36 s 


91 


+41°19'29'.'6 


2010/10/29.478 


17.8(U) 


2010/10/29.213 


Fell 


23 


M31N2010-10e 4 


00 ,i 42 m 57 s 


76 


+41°08'12'.'3 


2010/10/31.727 


18.1(R) 


2010/11/01.158 


He/N 


24 


M31N2010-1M 5 


00' , 42 m 43 s 


58 


+41°12'42'.'6 


2010/10/11.10 




2010/10/11.110 




25 


M31N2010-lla 


00' 1 42 m 31 s 


60 


+41°09'51'.'5 


2010/11/23.05 


18.2(Swift uvwl) 


2010/12/01.250 


He/N 


26 


M31N2010-12a 6 


00 ,, 40 m 24 s 


37 


+41°04'03'.'5 


2010/10/01.15 


16.7(U) 


2010/12/01.270 


Fell 


27 


M31N2010-12b 


00 ,I 42 m 31 s 


08 


+41°27'20'.'3 


2010/12/11.095 


16.7(U) 


2010/12/10.115 




28 


M31N2010-12c 3 


00 h 42 m 5 g S 


67 


+41°17'21'.'2 


2010/12/11.095 


17.2(U) 


2010/12/15.093 


He/N 


29 



Re fer ences. — (1) [Rodriguez-Gil ct al. 2009; Ovcha roTeFall [20091; (21 IQvcharov et al.ll2 009; Hcnze ct al. 2009b: Mcdv edev eFall 
2009; (3) Hcnze ct al. 2009a; Barsukova ct al. 2009b; (4) Fabrika et al. 2009b; (5) Nakano & Yusa 2009; Podigachoski ct al. 2009b; 
Hornoch & Kusnirak 2009; Fabrika ct al. 2009a; Hornoch ct al. 2009a; Nakano ct al. 2009; Nakano & Itagaki 2009; Yamaoka ct al. 
2009; Martignoni et al. 2009; Di Millc ct al. 2009; Barsukova ct al. 2009a; Hcnze et al. 2009b; (6) Nishivama k, Kabashima 2009a; 
Bur witz et al.l 120091; IHornoch et al.l I2009CT (7) iNishivama et alj l2009a ; Hcnze et alTf20b 9d; Kasliwal 2009; Hcn ze et aU 12 009c: 
Kasliwal et al. 2009; (8) Nishivama & Kabashima 2009b; Podigachoski ct al. 2009d; Hornoch ct al. 2009b; Podigachoski ct al. 2009c; (9) 
Green 2009; Pictsch ct al. 2009b; Cagas 2009; Kusnirak ct al. 2009; Hornoch & Pcicha 2009; Pictsch ct al. 2009a; (10) Nishivama ct al. 
2009bj; (11) iBurwitz et all 1201 0c; Horiio chet all I2010M4 JGreerJ I2010el: IPietsch fc HenzeTTaoi Oa: Hor noch et all 120104 IBurwitz et all 
2010d; Pietsch & Henzc 2010b; (12) Hornoch ct al. 2010c la i (TlHenze et alll2010d; I Pietsch et alH2010il; (13 ) IHornoch et al.ll2010bllalli|; 
(14) Burwitz ct al. 2010a; Pictsch ct al. 2010d; Hornoch 2010; Hornoch ct al. 2010a c; Burwitz ct al. 2010b; (15) Pictsch ct al. 2010c; 
Hornoch ct al. 2010i; Pictsch ct al. 2010f; Hcnze ct al. 2010c; Barsukova ct al. 2010; Pictsch ct al. 2010i; (16) Hornoch & Zaschc 2010; 
Henze et all I2010al; IHornoch et al l I2010allil; IHenze et all I2010blid; IPietsch et all 1201011; (17) lYusal I2010a|; IShafter et al. 2010d; (18) 
Yusa 2010b; Pictsch ct al. 2010k; Shafter ct al. 2010f h; (19) Romadan ct al. 2010; Shafter ct al. 2010c; Krushinski ct al. 2010; (20) 
Yusa 2010c; Shafter ct al. 2010h; Pictsch ct al. 2010i; (21) Corral-Santana ct al. 2010; Shafter ct al. 2010i; (22) Nakano & Yusa 2010; 
Shafter ct al. 2010g; (23) Nishivama & Kabashima 2010; Sun ct al. 2010; Hornoch ct al. 2010d; Homochova & Wolf 2010; Shafter ct al. 
2010b; (24) Homochova ct al. 2010; Pictsch ct al. 2010b; Shafter ct al. 2010a; Pictsch ct al. 2010g (25) Pictsch ct al. 2010i; Hcnze ct al. 
2010d; Shafter ct al. 2010c; (26) Pictsch ct al. 2010h; Cao ct al. 2010; Fabrika ct al. 2010; (27) Zheng ct al. 2010; (28) Green 2010a; 
Koishikawa 2010; Pietsch ct al. 2010a; Green 2010b c: Yusa 2010d: Koishikawa ct al. 2010; Pietsch ct al. 2010c; (29) Pricto ct al. 2010a; 
Yusall2010d ; IG7ee n 2010d f ; IPrieto et al.ll2010bl ; 

1 The bracket following the magnitude indicates the observation band and U represents unfiltered observation. 

2 RN; 1997-llk=2001-12b=2009-llb. 

3 According to CBET 2610, 2010-Ola and 2010-12c are two different objects within 0".8. 

4 RN; 1963-09c=1968-09a=2001-07b=2010-10e. 

5 The nova 2010-10f is probably inside the Globular Cluster Bol 126. 

6 Another recurrent burst of PT And. 
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TABLE 3 
Novae observed with GALEX 



Name 


Epochs 


UT start 


UT end 


M31N 2009-09a 


63 


2009/08/20.429 


2009/10/15.707 


M31N 2009-10a 


13 


2009/10/02.629 


2009/10/21.457 


M31N 2009-10b 


22 


2009/10/08.587 


2009/11/15.725 


M31N 2009-llb 


63 


2009/08/20.429 


2009/11/15.725 


M31N 2010-06a 


13 


2010/08/12.894 


2010/08/21.998 


M31N 2010-06b 


13 


2010/08/12.894 


2010/08/21.998 


M31N 2010-06d 


13 


2010/08/12.894 


2010/08/21.998 


M31N 2010-07a 


13 


2010/08/12.894 


2010/08/21.998 



TABLE 4 

X-ray and Spectroscopic observation epochs of M31 novae 





Observation UT 




Novae 


Spectrum X-Ray 


Reference 



M31N 2009-08e 
M31N 2009-09a 
M31N2009-10a 
M31N2009-10b 



M31N2009-lla 
M31N 2009-llb 

M31N2009-lld 
M31N2009-llc 
M31N 2010-Ola 

M31N 2010-06a 
M31N2010-06b 
M31N2010-06d 
M31N 2010-07a 
M31N2010-09a 
M31N2010-09b 

M31N 2010-09c 
M31N2010-10a 
M31N2010-10b 



M31N2010-10c 
M31N2010-10d 
M31N2010-10e 



2009/09/18.014 
2009/09/18.18 
2009/10/09.91 
2009/10/09.820 

2009/10/15.5 
2009/10/19.936 
2009/11/25.32 
2009/11/08.35 
2009/11/09.82 
2009/11/21.31 
2009/11/22.19 
2010/01/13.08 
2010/01/15.09 
2010/07/19 
2010/07/19 
2010/08/08.86 

2010/07/19 
2010/09/07.23 
2010/10/01.39 
2010/10/06.40 
2010/09/14.23 
2010/10/06.38 
2010/10/28.31 
2010/11/03.31 
2010/11/11.28 
2010/10/17.35 
2010/10/30.31 
2010/11/04.32 



M31N2010-10f 2010/11/13.285 



M31N2010-lla 
M31N2010-12a 
M31N2010-12c 



2010/12/05.74 
2010/12/30.24 
2010/12/30.25 



2010/11/15.16 
2010/11/17.09 

2010/11/03.04 
2010/11/06.20 
2010/11/07.06 



Mcd vedev et al. 2009 
Barsukova ct al. 2009b 
Fabrika et al. 2009b 
Fabrika et al. 2009a 
Yamaoka et al. 2009 
Barsukova ct al. 2009a 
Hornoch et al. 2009c 

Kasliwal 2009 
Kasliwal et al. 2009 
Hornoch et al. 2009b 
Hornoch fc Peicha 2009 
Hornoch et al. 2010h 
Hornoch et al. *2010g 
Hornoch et al. '2010c 
Hornoch et al. °2010c 
Barsukova ct al. 2010 
Hornoch et al. 2010e 
Sha tter et al. 2010d 
Shafter et al. 2010f 
Shafter et al. 2010h 
Shafter et al. 2010c 
Shafter et al. "OlOh 
Shafter et al. 2010i 
Shafter et al. 2010i 
Shafter et al. 2010i 
Shafter et al. 2010e 
Shafter et al. ~010b 
Shafter et al. "2010a 



Pietscli^^a 
Pictsch ct a 



Shafter et a^ 



2010k 
2010b 



2010e 



Pictsch ct al. 2010i 
Pietsch et al. 2010i 
Pietsch et al. 2010i 
Fabrika et al. 2010 
this paper 
this paper 



All three methods show consistent magnitude measure- 
ments. Hence we adopt the first method. 

2.2.3. Photometric and Spectroscopic Follow-up 

Photometric follow-up of M31N 201 0-lla was triggered 
for the Palomar 60-inch telescope (iCenko et al.l 120061) 
on 2010 December 4. On 2010 December 30, we also 
used the Low R esolution Imaging Spectrometer (LRIS; 
lOke et all 119951 ) mounted on the Keck I telescope to 
obtain spectra of M31N 2010-12a and 2010-12c. These 
spectra were obtained with a configuration of the 560 
dichroic, the 400/3400 grating in the blue side and the 
400/8500 grating in the red side with a central wave- 
length of 7800 A. 



2.3. GALEX 

GALEX is a wide-field imaging UV space telescope 
that was originally launched with near-UV (NUV and 
far-UV (FUV) detectors. Unfortunately, the FUV detec- 
tor failed before the start of the PTF survey of M31 The 
NV detector effective wavelength is 2316 A and the band 
width is 1060 A (jMartin et all 120051: iMorrissev et al.l 
120071). This channel has a FoV of 1.2 square degree 
(|Morrissev et al.ll2007t) . 

2.3.1. GALEX Observation Summary 

The visibility of M31 for GALEX is defined by the Sun, 
Moon and Earth-limb constraints that limit the obser- 
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Fig. 2.— The temporal coverage of both PTF and GALEX. A 
plus sign indicates an epoch observed by both PTF and GALEX 
while a cross indicates epochs observed by neither. PTF-only ob- 
servations are denoted with circles. GALEX-only observations are 
represented by diamonds. Statistically, during the first observing 
season, PTF observed M31 for 20% of 156 total nights, GALEX ob- 
served for 21% nights and the two overlapped on 6% nights. In the 
second season, PTF observed for 75% of 169 nights and GALEX 
overlapped with PTF on nine nights. 

vations more than for ground-based observatories. The 
result is that GALEX was able to intensively monitor 
M31 in several campaigns that overlap with the PTF 
campaigns of 2009 and 2010 (2009 August 08 to 2009 
September 04; 2009 October 02 to 2009 October 21; 2009 
November 15 to 2009 November 22; 2010 August 12 to 
2010 August 21). Figure [5] illustrates GALEX temporal 
coverage and its overlap with PTF. A field centered in 
M31 is monitored daily, while ten outskirt fields in M31 
are observed every other day. So we have fifty-one epochs 
of the center field while about 22—27 epochs of the out- 
skirt fields. The exposures were all taken during single 
orbit pointings and thus have durations close to 1500 
seconds with a range from 740 seconds to 1703 seconds. 
Table Q] summarizes parameters of observations. 

2.3.2. Photometry 

Photometry was done with an aperture with a radius of 
six arcseconds at each nova position, and was calibrated 
using the stan dard GALEX zero-poin ts to the AB mag- 
nitude system (Morriss ev et al.l l2007). Eight novae were 
detected; see Tabled 

2.4. Data from Other Sources 

For completeness in optical and UV bands, we incor- 
porate photometric measurements of novae from ATel 
and CBET reports (see the reference column in Table 
[2|). This includes optical photometry in R, B, V, r', i' 
and Ha bands as well as unfilte r ed dat a. We also include 
data published in lShafter et al.l (|2011af ) where several no- 
vae in 2009 are well-sampled. The UV data in ATel re - 
ports are mainly obtained by Swift (|Gehrels et al.l l2004). 
For consistency of data between GALEX and Swift, all 
Swift data are calib rated onto the AB magnitude system 
(jSiegel et al.l 120101 ). When data from different sources 
are used, due to difference of filters, data may have sys- 
tematical errors among different projects. 

In addition, we include X-ray and spectroscopic infor- 
mation of novae in M31 from Swift from either ATel re- 
ports or our own analysis (see Table SJ. 



2.5. Missed Novae 

During the two seasons presented here, we did not find 
three announced novae in the PTF archive: M31N2009- 
10c, M31N2009-llc and M31N 2010-07b. 

M31N2009-10c took place at a = 00 /l 42 m 45 s 8, 5 = 
+41 15'57'.' (J2000), less than thirty arcseconds from the 
center of M31. This nova r eached maximum light, i.e ., 
R = 17.2, on Oct 9, 2009 (|Podigachoski et alJl2009ah . 
PTF observed M31 on Oct 10, 2009. In the vicinity 
of the nova position, if we fit a constant background, 
the mean count of ~ 50000 with a standard deviation of 
a ~ 6000. Given that the seeing then was about three 
arcseconds and that the zeropoint is 27.6, the 3-c limiting 
magnitude is 15.2. If the background is fit with a plane, 
the standard deviation is roughly 1500 and the limiting 
magnitude is 16.8. Thus we did not detect this nova in 
the images. In the subtraction image, the circle of thirty 
arcseconds in radius in the center of M31 is excluded in 
the subtraction algorithm where the algorithm does not 
work well. As a result, we missed this nova in our PTF 
data. 

There are gaps of forty arcseconds between adjacent 
CCDs. M31N2009-llc unfortunately fell in the gap be- 
tween two adjacent CCDs. 

M31N 2010-07b had a peak luminosity o i R = 20.7 mag 
on 20 10/06/04.135 and slowly decayed (Hornoch et al. 
I2010kft . which is below the single- image detection thresh- 
old of PTF. The second observing season of PTF started 
at the end of June (see Figure [T]), so we do not have 
enough images to co-add for detection of this nova. 

Admittedly, we did not carry out a real-time search 
of novae in M31 in the first two seasons of PTF. Our 
nova sample inherits the sensitivity and incompleteness 
of those projects that discovered these novae. 

2.6. Special Novae 

2.6.1. M312010-09c 

This nova has a faint and close neighbor with R= 
19.85 ± 0.2 within two arcseconds. In SDSS Data Re- 
lease 8, we found two nearby sourc es within two arc - 
seconds. After color transformation (jJordi et al J 12006). 
the neighbors have R= 20.53 and R= 20.33, giving a to- 
tal magnitude of 19.67 that is consistent with our mea- 
surements. The uncertainty caused by the neighbors 
has been included in the nova light curve. The spec- 
troscopic confirmation in ([Shatter et alJl2010cft identified 
the M31N 2010-09c as a real nova of Fe II type. 

2.6.2. M312010-10} 

This nova is spati ally associated wit h the M31 globu- 
lar cluster Bol 126 (jWirth et alj 1 1985ft . In Figure H we 
clearly see a brightening of the globular cluster by about 
one magnitude on around 2010 October 10 (Julian Date: 
2455480). The luminosity of the globular cluster before 
the nova is 16.8 ± 0.2. The light curve of the nova is 
obtained by subtracting the globular cluster brightness. 
One possible reason for the large uncertainty in the pho- 
tometry measurements is that the nova is only a few 
pixels away from the edge of the CCD. 

2.6.3. M31N2010-01a and M31N 2010-12c 

These two novae occurred within one arcsecond of each 
other, a separation which is less than the typical spatial 
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resolution of PTF (~l'/0). A highly accurate astrometric 
solution show s that they are indeed two separate objects 
(|Greenll2010d) . 

We obtained a spectrum of 2010-12c with LRIS on 2010 
December 30 (Figure 02 bottom panel). The spectrum 
shows a weak continuum superposed with many emission 
lines, illustrating that the nova has already entered the 
nebular phase. We identified Ha, H/3, H7, H5, 01 8446 
and several He and N lines. H5 has a slight P Cygni 
profile. The emission lines suggest the nova is of the 
He/N type. After fitting the emission lines with a Gaus- 
sian profile, we find that Ha is centered at 6556 A, H/3 
at 4858 A, H<5 at 4338 A, and OI 8446 at 8437 A, sug- 
gesting a blueshifted velocity in the range from —200 to 
—300 km s -1 . On the other hand, the systemic velocity 
of M31 is —300 km s . The nova lies abo ut 2' northeast 
of the center of M31. Th e HI 21 cm line (jChemin et al.l 
120091 : ICorbelli et al.ll20l"ol ) suggests a rotation velocity of 
100 km s _1 receding from us. However, the nova is ap- 
parently close to the center of the galaxy, and therefore 
may not share the motion of the disk measured by the 
HI 21cm line. Another measurem ent of a nearby plan- 
etary nebula (|Merrett et al.ll200ll at a = 00 ?l 42 m 57 s .4 
and 5 = +41 17'26'.' (J2000) suggests a radial velocity of 
—300 km s _1 along the line of sight. Consequently, we 
conclude that 20 10- 12c is a real nova in M31. 

Additionally, these two novae have two more histori- 
cal neighbo ring novae within 10": nova 30 (discovered in 
1986) from ICiardullo et al.1 (Il98l and M31N 2009-08a. 
Generally speaking, we observed more novae in the cen- 
tral region of M31 than in outskirts. This might result 
from the high stellar density in the center of the galaxy. 

3. LIGHT CURVES OF CLASSICAL NOVAE 

Except three recurrent events, light curves of twenty- 
six classical novae are presented in Figure 03 to 03 An 
electronic version of the PTF and GALEX data can be 
found onlinj^l in a machine- readable format. An exam- 
ple of M31N2010-12a is illustrated in Table 03 

The nova light curves are presented in groups based 
on their features described below. Figures 03 presents no- 
vae with smooth decline. Figure 03 illustrates novae with 
jittery decline. Novae with well-sampled rise stages ap- 
pear in Figure 00 Figure 03 highlights two novae recorded 
during the PTF lOmin-cadence observations where their 
rise histories were recorded in great detail. Figure OS 
shows the light curves of the only globular cluster nova 
M31N2010-10f and its host Bol 126. All other under- 
sampled light curves are collected in Figure 03 

Table 03 summarizes the photometric and morphologi- 
cal features of the novae. 

3.1. Morphological classification 

The shape of each light curve can be characterized by 
several parameters: the rise rate, the rise time to the 
maximum magnitude, the rise behavior, the maximum 
magnitude, the decline rate that is represented by t n (n = 
1, 2, 3, . . .) defined as the time interval during which 
a nova decays from its maximum by n magnitudes, and 
the decline behavior. Since PTF samples the M31 nova 
luminosity distribution down to ~ —4.5 (see §1) and 

13 http:/ /www. astro. caltech.edu/~ycao/m31novae/index. html 



TABLE 5 
Photometric measurements of 
M31N2010-12A 



Julian Date Magnitude Magnitude Error 



2455530 


.702 


20. 


.190 


99.000 


2455531 


.770 


16. 


.763 


0.030 


2455532 


.780 


15. 


.581 


0.023 


2455532 


.825 


15. 


.561 


0.025 


2455535 


.681 


16. 


.013 


0.024 


2455537 


.578 


16. 


.325 


0.024 


2455537 


.622 


16. 


.346 


0.025 


2455538 


.577 


16. 


.526 


0.024 


2455538 


.621 


16. 


.496 


0.025 


2455539 


.577 


16. 


.636 


0.026 


2455539 


.620 


16. 


.638 


0.025 


2455540 


.615 


16. 


.721 


0.027 


2455540 


.659 


16. 


.776 


0.026 


2455541 


.613 


16. 


.817 


0.027 


2455541 


.656 


16. 


.899 


0.031 


2455544 


.593 


17. 


.243 


0.029 


2455544 


.639 


17. 


.175 


0.032 


2455545 


.593 


17. 


.283 


0.038 


2455545 


.637 


17. 


.303 


0.034 



1 99.000 in Magnitude Error indicates that a 3(T 
upper limit is given in Magnitude. 

the peak absolute magnitude of a nova can range from 
—5 to —9, we use ti to indicate the decline rate. In 
our measurements, ti is estimated by linear interpolation 
between adjacent measurements. Other parameters are 
defined as follows: the rise rate is the mean rate from 
the first detection to the first local maximum; the rise 
time is the interval from first detection in PTF to its 
global maximum. All measurements are done with PTF 
R-band data. 

For novae inside our Galaxy iStrope et all ((2010) clas- 
sified the diverse light curves morphologically according 
to their evolution of decline. With a median coverage 
extending to 8.0 mag below peak, the authors grouped 
nova light curves into seven classes. In their sample of 
93 nova light curves, 38% decline smoothly, 21% have 
plateau phase, 18% show dust dips, 1% rebrighten with 
a cusp-shaped secondary maxima, 4% oscillate quasi- 
sinusoidally on smooth decline, 2% are flat-topped and 
16% have irregular jitters or flared during decline. How- 
ever, PTF is unable to follow M31 novae to such late 
times in the light curve evolution. We, instead, roughly 
classify our M31 novae into two major classes: S-class for 
smooth light curves and J-class for jittery light curves 
or for light curves with flares superposed on the decline. 
Statistically in our nova sample, 60% show S-class decline 
and 40% belong to the J-class. Given that our sample 
is quite small, this result is roughly consistent with the 
Galactic population. We do not see any flat-topped nova. 

Similarly, due to the high cadences, PTF sampled the 
rising stages of several novae very well. This enables 
us to employ the same S-J classification to describe the 
rising parts of the light curves (see Table 03) ■ We find 
that the S-class novae rise very fast (e.g., M31N 2010-09b 
rose within two days) while the J-class novae rise more 
slowly (e.g., M31N 2010-10b took thirty days to reach its 
maximum). But we do not see any relation between the 
rise time and peak magnitudes. Luminous M31N2010- 
09b and subluminous M31N 2009-lla rose rapidly while 
luminous M31N 2010-07a and subluminous M31N2009- 
09a both rose very slowly. 
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Fig. 3.— Spectra of M31N2010-12a (top) and M31N2010-12c (bottom), obtained with LRIS on Keck-I telescope on December 30th, 
2010. 
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Fig. 4. — The light curve of the globular cluster Bol 126 is shown 
in the top panel. The subtracted light curve of M31N 2010-10f is 
presented in the bottom panel. 



During a very-high-cadence (tens of minutes) ex- 
periment in PTF, we observed the rising stages of 
M31N 2010-09a and M31N 2010-09c in unprecedented de- 
tail, as shown in the bottom panels of Figure HI The rise 
of both novae is smooth over a timescale of tens of min- 
utes. 

3.2. Super- Eddington Phase 

Some novae around maximu m exhib it a super- 

Eddington phase. Theorists (e.g.. IShavivl I200TT ) model 
this phase with a porous structure in the nova wind that 
reduces t he effective opacity and raises the Eddington lu- 
minosity. iKato fc Hachisul (|2005|) produced a light-curve 
model based on this idea that successfully matched the 
super-Eddington phase of Nova V1974 Cyg, and demon- 
strated the importance of the relationship between the 



optical and UV light curves. They isolated a portion of 
the UV continuum near 1455 A that had a duration that 
was well correlated with the decline rate of the nova in 
the optical. This continuum UV light curve combined 
with the model of how the UV light is reprocessed into 
the optical light curve allows an estimate of the precur- 
sor WD mass and distance to be made, based on certain 
assump tions about the chemi cal composition of the en- 
velope (|Kato &: Hachisull2007l) . 

One might assume that the clumpy wind responsible 
for the lower opacity and higher Eddington luminosi- 
ties of the brightest novae might produce other observ- 
able signatures. Small-scale structure in the rising light 
curve could be observed as the clumps are non-uniformly 
heated and accelerated by the ongoing eruption and as 
radiation is diffused through them. We have very high- 
cadence observations of the rising phases of two novae 
(2010-09a and 2010-09c), both of which have a peak 
brightness fainter than -7.5 and thus are probably not 
super-Eddington novae. Nonetheless, while 2010-09c ex- 
hibits a very smooth and steady rise and a smooth de- 
cline, consistent with a uniform nova wind, 2010-09a has 
considerable structure in its light curve with multiple 
peaks both near and after the maximum light. Our ob- 
servations of 2009-09a, 2010-07a, 2009-09a and 2009-09c 
also show several peaks around maxima. The multiple 
peaks could arise when different winds are blown out 
from the center successively during an outburst. Each 
peak is formed when the pseudo-photosphere meets the 
shock wave at the head of the corresponding wind. We 
conclude that structures in the nova wind and thus in 
the light curve can arise regardless of the luminosity of 
the outburst. 

4. NOVAE IN GLOBULAR CLUSTERS 

M31N2010-10f was spatially associated with the glob- 
ular cluster (GC) Bol 126 to within one arcsecond. A 
supersoft X-ray source ( SSS) w as also reported in the 
vicinity by iHenze et ah (2010d). This event is similar 
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Fig. 5. — Light curves of nine smoothly declining novae. Different symbol shapes mean data from different sources. PTF and GALEX 
data are s hown in filled circles and rectangles, respectively. Data collected from ATel and CBET reports are in empty triangles and data 
read from Shatter ct al. (2011b) are in empty rectangles. Data in different filters are illustrated by colors. R band data are in red; B band 
data in blue; V band data in green; r' band data in orange; i' band data in Chartreuse; Ho narrow band data in cyan; NUV band data in 
magenta; and unfiltered data in black. In the top axis, an "X" indicates an X-ray observation at that time and an "S" indicates a spectrum 
of the nova taken at that time. 



to M31N2007-06b ((Shatter fe Quimbvl 12001 in Bol 111 
which was confirmed spectroscopically during the out- 
burst. A SSS count erpart was also found for this object 
(jHenze et al.ll2009eD . They are the only two events in 
M31 that are associated with GCs. 

Severa l previous searches for GC novae obtained nu ll 
results (|Ciardullo et al.1 Il990bt iTomanev et~aT1 11992( 1 . 
We examined our PTF image database at the positions 
of all confirmed M31 GCs from the Revis ed Bologna Cat- 
alog (V.4, Dec 2009. iGalleti et al.|[25ol . Magnitudes of 
these GCs range from R=13mag to R=20mag. We do 
not find any other GC novae candidates. Given an effec- 
tive survey time of M31 in PTF approximating one year, 
we obtain a GC nova rate of about 1 yr _1 . 

Is the nova rate enhanced i n GCs? The tota l stellar 
mass of M31 is ~ 6 x 1O 1O M (jTamm et al.ll2007l) . Using 
the highest pu b lished nova rate for M31 of 65 yr _1 from 
iDarnlev et al.l (|2006l ). we estimate a maximum global 
specific nova rate of ~ 1/(10 9 M Q ) yr _1 in M31. Given 
a typical GC mass of 10 6 Af© and about seven hundred 
confirmed M31 GCs, the nova rate in the M31 GC sys- 



tem is ~ 1/(10 9 M Q ) yr" 1 . We conclude that the M31 
GC specific nova rate is not significantly higher than the 
M31 galaxy specific nova rate. 

LMXBs - mass transferring binaries with a neu- 
tron star receiving matter from a companion - can be 
viewed, in some ways, as surrogates for CNe and RNe. 
Since the early days of X-ray astronomy it has been 
known that at least in our Galaxy that the brightest 
LMXBs (per unit mass) are nearly two orders of mag- 
nitude larger than that for the disk of the Milky Way. 
This trend appears to be even more acute for M31 
(|Trudolvubov fc Pr icdhorskv 2004). In part this may be 
because M31 possesses a more extensive cluster system 
than our own Galaxy. Thus, naively, we would expect 
far more CNe and RNe in the globular cluster system of 
M31. Over the period of investigation (2009-2010), a to- 
tal of three RNe were found in the disk/bulge of M31. If 
we assume the mass fraction of globular clusters is ~ 1%, 
we would have expected to see a similar number in the 
cluster system. Within small number systems, the detec- 
tion of RN associated with the cluster would be consis- 
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Fig. 6. — Light curves of novae with jittering decay. Symbol designation is the same as that in Figure [5] 
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Fig. 7. — Light curves whose rise stages are well-sampled. Symbol designation is the same as that in Figure [5] 
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Fig. 8. — M31N 2010-09a and M31N 2010-09c were observed during the PTF lOmin-cadence experiments. In the top two panels, we show 
both full light curves. In the bottom two panels, rising phases are highlighted. Symbol designation is the same as that in Figure(5] 
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tent with this expectation. Separately, given that seven 
hundred CNe are known, over 2009 and 2010 we should 
have seen many tens of CNe in the globular clusters, 
but only one was found. Indeed, over the last hundred 
years only one CN that has been associated with a glob- 
ular cluster of M31 has b een spectroscopically confirmed 
(jShafter fe Quimbvl[2007h . In summary we cannot arrive 
at a sensible conclusion with the presently available data. 

5. RECURRENT NOVAE 

According to the catalog maintained by lPietschl (2010), 
26 RNe in M31 are known. During the two seasons 
we report on, three recurrent events were reported: 
M31N2009-llb, M31N2010-10e and M31N 2010-12a. 
Figure [TU1 presents their light curves. Besides these, we 
examined the PTF archive at the posit ions of more th an 
800 old novae in M31 from the list of lPietschl (|2010l) in 
order to find more recurrent events. None were found. 
Therefore, given an effective observing time of roughly 
one year in PTF, the recurrent nova rate in M31 is 
~ 3 novae per year. 

We examined M31N 2010-12a in some detail, as 
another outburst of PT And. This object has 
recorded outbursts i n 195 7, 1983, 1986, 1988 and 1998. 
lAlksnis fe Zharoval (|2000D collected its previous out- 
bursts and found that all light curves can be fit by one 
template. The average decline rate of previous outbursts 
is 0.09 to 0.11 mag per day. The light curve of this most 
recent outburst (M31N 2010-12a) can also be fit by the 
template. Figure QT] shows the similarity between light 
curves of the recent and 1988 outbursts. The decline 
rate of the recent outburst is around 0.10 mag per day 
consistent with previous ones. Th is invariance p roperty 
is also observed in Galactic RNe (Schacfcr 2010). Such 
characteristics of RNe suggests that RN outbursts only 
depend on system parameters like the WD mass and its 
composition. 

Because of its luminous and frequent outbursts, PT 
And was assumed to be a dw arf nova inside the Milky 
Way $Mk snis fe Zharoval I2000D . However, no direct ev- 
idence of this object's nature was available until we ob- 
tained a spectrum of its most recent outburst with LRIS 
in Keck-I telescope. This spectrum (see the top panel 
of Figure [3]) shows a very blue continuum superposed 
with prominent Ha and H/3 lines as well as several weak 
Fe lines. Fitting the red part of the spectrum with a 
Planck function yields a temperature of 7 ~ 8 x 10 3 K . 
We fit the Ha and H/3 lines with Gaussian profiles and 
obtained the following observed wavelengths: Ha is cen- 
tered at A = 6558 Aand H/3 at A = 4857 A. These corre- 
spond to blue-shifts with radial velocities of ~ —200 km 
s^ 1 and ~ —300 km s _1 . The systemic velocity of M31 
is ~ —300 km s _1 . We assume that this object is not 
located in the disk of M31, thus we do not take the rota- 
tion velocity of M31 disk into account. We conclude that 
the blue-shifts are roughly consistent with the motion of 
M31. Our spectrum evidently supports the conclusion 
that PT And resides in M31 instead of in the Milky Way. 
This can also be tested by observing its proper motion. 
It is of only galactic sources that we may observe proper 
motion. In addition, the existence of weak Fe lines sug- 
gests the nova is of Fe II type. 

The association of PT And with M31 leads to a peak 
absolute magnitude of about -9, which makes it the 



brightest RN ever known. 

The decline rates of RNe vary over a large range in 
M31. M31N2009-llb has t 2 > 37 days and M31N2010- 
lOe has t 2 > 5 days. Consider that 2010-12a has a decline 
rate of ~ 0.1 mag/day, so its t 2 is possibly around twenty 
days. The t 2 valu es of Galactic R Ne also range from one 
day to fifty days (ISchaeferll2010h . 

6. UV LIGHT CURVES OF NOVAE 

Theoretically, the UV emission after the optical max- 
imum light is characterized by a delayed UV peak af- 
ter the optical peak, because most UV photons are ab- 
sorbed and re-emi t ted in t he optical in the outer envelope 
(iKato fe Hachisul Il994t [G onzalcz-Ries tra fe Krautterl 
119981: ICassatella et all I2002l ~ However in our sample, 
we see some deviation from this general picture: 2009- 
10b has a peak in the GALEX NUV over two days prior 
to its peak in the optical (see Figure [T2")) . The recurrent 
nova 2009-llb was seen in the GALEX NUV about forty 
days before it was visible in optical bands (although this 
could be due to cadence and weather, se e Figur e fTUl) . 
The n ova 2010-06d was detected in Swift (Pictsc h et al.l 
20101) about six days before optical detection (see Fig- 
ure |9]). Swift also detected 2010-lla several days befor e 
optical detection (jPietsch et al.ll2010bJ: iCao et alJfeOlOl ; 
also Figure [9]) . For 2010-07a, the optical and UV max- 
ima were nearly simultaneou s as measured by PTF and 
by Swift (jHenze et al.ll2010bl see Figure (H). 

Admittedly the detection in the optical and UV in- 
volves many factors like survey sensitivities, temporal 
coverage and observing weather. In a few cases, how- 
ever, it is clear that the UV peak coincides with or pre- 
cedes the optical peak. These cases demonstrate that not 
all UV photons are re-processed and that the absorption 
in the UV will vary considerably depending on factors 
such as physical geometry or chemical abundances in the 
envelope. These variations will produce a range of opti- 
cal/UV peak offsets and luminosity ratios. We also point 
out that multi-band optical data show that three of these 
UV nov ae have red colors: 2 009-10b has R — V ~ 1.4 at 
decline dShafter et al.ll2011bf): 2 009-llb has R - V ~ 0.7 
at decline (jShafter et alJl2011bt ); 2010-lla has g-i ~ 0.2 
around its optical peak from our photometric follow-up. 

6.1. M31N2009-10b and M31N 2010-lla 

The novae 2009-10b and 2010-lla are two peculiar 
cases in that their UV peaks lead the optical maxima 
by a few days. Moreover, at the early decay stage of UV 
emission, the optical colors of both novae look quite red. 

Of the two, 2009-10b is brighter, has better light curve 
coverage in PTF, GALEX and other bands, and is thus 
a better case for further examination. Figure [12] zooms 
in on the early multi-band light curve. Table |8] lists the 
rise rate, maximum magnitudes, dates and decline rates 
by t 2 and t 3 in the NUV, B, V and R bands. The NUV 
data are from GALEX. Part of R band data are from 
PTF. B, V and part of R band data are published in 
IShafter et all (|2011bf ). We can see clearly that the NUV 
peak is two-to-three days earlier than the B, V and R 
peaks. The rise rates in B, V and R band are consistent 
while NUV emission rises swiftly. The decline rates of 
NUV, B and V are similar while R band emission decays 
slower. 
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Fig. 10. — Light curves of 3 RNe detected in PTF. Symbol designation is the same as that in Figure [5] 



TABLE 6 

Properties of PTF well-sampled nova light curves 
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The maximum magnitude is from ATel reports because PTF missed the peak magnitude. 
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TABLE 7 

P60 PHOTOMETRIC FOLLOW-UP OF M31N 2010-llA 
JD-2455500(days) Bands Magnitude Magnitude Error 



45 50 55 60 65 70 75 80 85 90 

JD-T 

Fig. 11. — A comparison of the light curves of outbursts of PT 
And in 1988 and 2010 is shown in the figure, where T = 2455466 for 
the 2010 outburst and T = 2455100 for the 1988 one. Though with 
different filters, the light curve of the recent outburst (M31N 2010- 
12a) resembles that of the 1988 outburst. 



2010-lla is possibly the faintest nova in our sample 
with a peak observed magnitude of R ~ 18.8 (see Fig- 
ure 0). It was first detected by SWIFT (jPietsch et al.l 
2010 h) and the n we reported its optical counterpart 
(|Cao et aLll2010h . The strong UV emission peaked sev- 
eral days prior to its optical maximum. Near its optical 
maximum, we carried out a photometric follow-up in i', 
r', g' and B bands with the Palomar 60-inch telescope 
and observed that this nova had red optical colors, as 
shown in Table [7] 

The UV observations and the color information sug- 
gest that UV and optical emission originates from dif- 
ferent radiative regions because they cannot be simply 
interpreted as the thermal radiation of a receding spher- 
ical photosphere which has a large optical depth in local 
thermal equilibrium. Because UV emission is only gen- 
erated deep inside the envelope, the result may suggest 
that the envelope has holes from which we can see deep 
inside. 

The models of lShavivl (l200l : IKato fc HachTsul (l200l 
2007) propose a reduced effective opacity of the nova en- 
velope to explain the super-Eddington phase, resulting 
from a porous structure of the envelope caused by its 
fluid instability. We may use a similar scenario of as- 
phericity of nova outbursts to qualitatively explain the 
red UV novae. UV photons are generated deep under 
the envelope at the surface of the WD. In the envelope, 
the bottom is heated to very high temperatures by ther- 
monuclear runaway, and thus instabilities like Rayleigh- 
Benard convection are likely to occur. This instability 
causes the mass distribution in the nova envelope to de- 
viate from a spherical distribution. Thus, along some 
lines of sight to the center of the outburst, the mass is 
less and so is the opacity. The UV photons then escape 
from these directions and are observed by us if we are ob- 
serving along one of these lines of sight. The red optical 
emission still comes from the reprocessing of the enve- 
lope in the other directions. The ratio between the UV 
and optical fluxes reflects the angular size of the cone. 
If the instability produces an extreme feature with very 
low mass and very low opacity along the line of sight, we 
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Fig. 12.— The multiband light curve of M31 N2009-10b focusing 
on its early evolution. We can clearly see that the UV peaks prior 
to the peaks in B, V and R bands by about two days. When the 
UV decays, the optical part is quite red. 

can observe the central UV photons directly. This may 
be what we are observing in 2009-10b and 2010-lla. An- 
other possible explanation could be very low abundances 
of chemical species that absorb in the UV, although this 
may not explain the red optical colors. 

The uniqueness of the UV-optical light curves of 2009- 
10b and 2010-lla suggest either their super-Eddington 
phases or their chemical abundances (o r both) are very 
differ ent from the majority of novae ([Cassatcll a et al.1 
2002). They may possibly represent a population that 
deviates from any empirical relation derived from large 
nova samples, such as MMRD relation. In fact, 2009- 
10b is a luminous outlier of the MMRD relation (see 
Figure IT3| . Unfortunately, due to its faintness, we were 
unable to measure the ti decline rate for the nova 2010- 
lla. 

7. MMRD REVISITED 

Novae exhibit a relation between t he maximum magni- 
tude and the rate of declin e (MMRD Idella Valle fc Liviol 
IT995t IKasliwal et al.ll27)Tol ). We can assess this relation 
with our nova sample in M31. 

One of main issues in evaluating the MMRD relation is 
estimating the dust extinction. A multitude of methods 
using color and Balmer decrement have been used in lit- 
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TABLE 8 

2009-IOb NUV and optical properties 



Bands NUV 1 B 2 V 2 R 3 

rise rate (mag/day) 1.2 0.4 0.4 0.7 

peak magnitude 16.68 15.24 15.00 14.84 

peak date (JD-2455110) 7.95 10.37 10.37 9.95 

t 2 (days) 3 5 5 8 

t 3 (days) 7 7 7 17 

1 GALEX o bservations; 

2 data from lShafter et al.l H201lbl) ; 

3 data from PTF and;Shaftcr ct all H2011bl) . 



erature to es t imate the extinction (|Darnlev et al.l [20061 : 
IKogurd[l96lt IKasliwal et alllMol) . The PTF nova ob- 
servations, however, provide no color information, being 
done only in the R band. Thus, we apply a mean fore- 
ground dust extinction. Th e foreground extinc tion of 
M31, E(B — V) = 0.062mag (|Schlegel et al.|[l998l) . leads 
to 



A g = 3.793 x E(B - V) = 0.234. 



(1) 



For a better compar ison with previous researches (e.g., 
IKasliwal et~al"1l2010l) . following Shatter et al.(2009), we 
use the colors of an A5V star (T = 8200-ff) to tra nsform 
R magnitudes to g mag nitudes (|Jordi et al.ll2006D . 

The MMRD relation is presented as the gray band 
in Figure [13] We overplot on this diagram our 
M31 nov ae (red points) and compa re them to the no- 
vae from Idella Valle &: Livid (|1995l) (gray points) and 
IKasliwal et al.l (|2010l) (gravstars). Without a local ex- 
tinction correcti on, most of our peak m agnitudes are un- 
derestimated. In Kasliwal ct al. (2010), the total extinc- 
tion correction for M31 novae is roughly A g — 0.24 except 
that two of them (20%) have a correction of A g > 1.0. 
It is possible that local extinction for most novae in our 
sample is small. 

In Figure 1131 we can see five significant outliers from 
the MMRD region in our sample. The brightest outlier 
is M31N2010-10b (see discussion in Section IBTTj) . The 
other luminous outlier is PT And (=M31N2010-12a). 
M31N2009-12a, M31N2010-06b and M31N2010-06c lie 
in the faint and fast-decline rate z one below the MMRD 
relation shown by the gray band. IKasliwal et al.l (120101) 
also f ound several novae in this zone. Kasliwa l et al.l 
(2010) proposed the possibility that the outliers may be 
RNe. However, none of the three faint outliers in our 
sample are known to be recurrent. The remaining novae 
reside roughly along the MMRD region. 

8. CONCLUSION 

In this paper, we tabulate the photometric measure- 
ments of twenty-nine classical novae in M31 from PTF 
data and present their optical light curves, eight of which 
also have joint UV-optical light curves from both PTF 
and GALEX. The main findings are as follows: 

1. We detected three recurrent novae. They show sim- 
i lar properties to recurrent novae in the Galaxy 
(jSchaeferl I2010D . We obtained a spectrum of 
the recent outburst of PT And (also known as 
M31N2010-12a). While previo usly identified as a 
dwarf nova in the Milky Way ([Alksnis fc Zharoval 
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Timescale "t2" [day] 

Fig. 13. — The empirical relation for novae between the max- 
imum magnitude and the rat e of decline (MMRD). Th e gray re- 
gion denotes the MMRD from ldella Valle fc Liviol l| 19951 ). while the 
small dark gray solid dots repr esents their nova sam ple. The six 
gray stars show the sample from Kasliwal ct al. (2010). Red circles 
and a red rectangle are the M31 novae from this paper. The red 
rectangle denotes the bright recurrent nova M31N 2010-12a. 



2000), we find the blueshifts in our spectrum to be 
consistent with the event being a recurrent nova in 
M31. The light curve of 20 10- 12a resembles those 
of previous outbursts, suggesting the same physical 
conditions in the recurrent outbursts. Our search 
at all other historical nova positions finds no addi- 
tional outburst candidates. From this we derive a 
recurrent nova rate of about three per year in M31. 

2. The rise behaviors of M31 classical novae are di- 
verse and we classify them morphologically into S 
(smooth) and J (jittering) classes. We find that 
the S-class novae rise significantly faster than the 
J-class novae. We see no dependence between rise- 
time and maximum magnitude. 

3. During the ten minute cadence PTF observations 
the rising light curve of two novae were recorded 
with unprecedented detail. These novae showed 
smooth and consistent rising light curves without 
large amplitude variability. 

4. The declining light curves are also divided morpho- 
logically into S and J classes. The relative popula- 
tion of different classes in M31 is roughly consistent 
with that in the Galaxy. 

5. Three quarters of the well-sampled novae obey the 
empirical MMRD relation. The remaining one 
quarter (five events) are co nsistent with the "fain t 
and fast" outliers found by IKasliwal et al.l ([2010). 

6. Some novae have UV detections prior to or at the 
same time as the first optical detection. Moreover, 
2009- 10b and 2010-1 la have UV peaks prior to op- 
tical maxima while both are quite red in the UV 
decay phases. These observations are inconsistent 
with theoretical predictions of delayed UV emis- 
sion in the iron curtain stage, possibly suggesting 
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extreme density or abundance variations for these 
novae. 



7. Apart from 2010-10f in Bol 126, we find no other 
novae in all the cataloged M31 GCs. We derive a 
GC nova rate in M31 of ~ 1 yr _1 . This nova rate 
is not enhanced relative to the rate in M31 outside 
of the GC system. 



We realize the novae reported here are only a subset of 
all the classical novae in M31 as t he number is smalle r 
than expected (e.g., 65 nova yr" - 1 IDarnlev et all [2006). 
We did not undertake a systematic search for classical 
novae in PTF images during the 2009 and 2010 seasons. 
Since the FoV and survey depth of PTF should be sensi- 
tive to all classical novae (weather-permitting), we have 
started a more complete search for classical novae in M31 
in the 2011 season. For 2012, we intend to run a real- 
time search for novae. Our photometric and spectro- 
scopic follow-up will enable us to better understand their 
properties and to estimate the local extinction correction. 

Besides optical bands, future UV surveys are war- 
ranted for nova studies. Some novae peak in the UV 
before the optical. Such novae could provide impor- 
tant insights into the development and structure of the 
nova wind. It is also important to know if novae like 
M31N2009-10b and M31N2010-lla are rare or common 
and if they are always outliers in the MMRD diagram. 

Ongoing and upcoming synoptic surveys (e.g., Palomar 
Transient Factory and Next Generation Transient Facil- 
ities) will continue to monitor the Andromeda Galaxy. 
The long time baseline will give us a larger sample of 
recurrent novae. As noted earlier, recurrent novae could 
be the progenitors of la supernovae. Thus the increased 
sample and timely follow-up will give us a comprehensive 
understanding of rate of recurrent novae-an important 
clue to their endpoint. 
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